Cracks in rubber under tension exceed the shear wave speed 
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The shear wave speed is an upper limit for the speed of cracks loaded in tension in linear elastic solids. We 
have discovered that in a non-linear material, cracks in tension (Mode I) exceed this sound speed, and travel 
in an intersonic range between shear and longitudinal wave speeds. The experiments are conducted in highly 
stretched sheets of rubber; intersonic cracks can be produced simply by popping a balloon. 
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Cracks advance by consuming the potential energy stored in 
the surrounding elastic fields and, thus, one expects that they 
should travel no faster than the speed of sound. While this 
intuition is indeed correct, solids support longitudinal, shear, 
and surface (Rayleigh) waves, with distinct speeds cl > c s > 
cr, respectively [1], and determining which of these is the 
upper limit for crack speeds is a subtle problem. 

In the case of linear elastic solids, Stroh was first to argue 
that cracks loaded in tension (Mode I) could not exceed the 
Rayleigh wave speed cr 12] ■ The mathematical analysis of 
Freund and others 1 3, 4] showed that the energy consumed by 
a crack diverges as its speed approaches cr, suggesting that 
cracks cannot travel faster than this speed. However, geologi- 
cal field measurements of fault motion |5], computer simula- 
tions 1 6, 7], and laboratory experiments |8] have established 
that cracks loaded in shear (Mode II) can exceed c s (and hence 
cr) and reach speeds close to ct,. Dynamic fracture theory 
was extended to include these "intersonic" cracks |4, 
However, even when revisited in light of the arguments that 
allowed intersonic cracks in shear-loaded samples, linear elas- 
tic fracture mechanics forbids intersonic cracks loaded in ten- 
sion 0. 

In the case of elastic materials subjected to large deforma- 
tions, it is natural to assume that the crack speed is bounded 
above by the shear wave speed. Recent numerical simula- 
tions have found mode I cracks running faster than this wave 
speed | llj. Here we present experimental results which show 
crack speeds faster than the shear wave speed occur in pop- 
ping rubber. 

Experiment. Our experiments are conducted with sheets of 
natural latex rubber 0.15 mm thick purchased from McMaster- 
Carr, stretched biaxially in a tensile testing machine, and 
punctured by pricking the sheet with a needle. The sam- 
ples are rectangular, 12.7 cm wide and 27.9 - 34.3 cm long, 
depending on the expected maximum extension, and their 
perimeter is divided into tabs 3.0 cm wide that serve as grip- 
ping points for the testing apparatus. A square grid is drawn 
on the sample to monitor the magnitude and homogeneity of 
extension as the sample is stretched. 

The apparatus consists of two fixed and two mobile linear 
guide rails (see Fig. [TJ a ))- The mobile rails move indepen- 
dently and orthogonally. Multiple sliding clamps that grip in- 
dividual tabs on the sample are mounted on each rail. The 
sheet was attached to these clamps and stretched to an exten- 



sion state {XxAy) such that X y > X x , where A; is the ratio of 
the stretched to original length along the i direction, and x and 
y are the long and short directions, respectively. As the rails 
are moved apart, the sheet expands and the clamps separate by 
sliding along the rail, yet remain equidistant due to the elastic 
coupling to each other through the sheet. 
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FIG. 1 : Apparatus for creating biaxial extensions in a thin rubber 
sheets, (a) Initial state of the sample loaded in the test apparatus, (b) 
Final state of the sample after extension. The dark rectangle is a rigid 
frame attached to the sample and the ' x ' marks the crack initiation. 



After reaching the desired extension state, the sheet is 
clamped between two steel rectangular frames (10.2 cm x 
66 cm)(Fig. ^b)). Then the extension state of the sheet is 
fixed and determined exclusively by these clamping frames. 
The frames fix the total energy available to the crack, and 
since the aspect ratio is small, the crack will consume equal 
amounts of energy per crack length 1 12]. 

A crack is initiated by pricking the sheet with a pin at the 
point marked ' x' in Fig.[Hb) and the velocity is determined 
photographically. Along the expected crack path we position 
an optical sensor that triggers a strobing flash when the crack 
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FIG. 2: Multiple exposure photograph of a crack in a rubber sample 
(X x = 2.2, X y = 3.2). The exposures, separated in time by 500 fis, 
are captured on a single frame with a CCD camera. 

parts the rubber blocking the beam path of the sensor. A CCD 
camera captures multiple images from each flash pulse on a 
single frame producing an image as shown in Fig. [2] As re- 
ported in 1 13], the crack can run straight or wavy depending 
on the extension state; irrespective of the crack path, the ve- 
locity of the crack, v c , is calculated by the distance between 
successive crack tip positions and the time interval between 
flashes. Cracks are wavy for \ x > 2.2. These measurements 
are reproducible within 5% from trial to trial. 

Longitudinal wave speed. We also measure the longitudinal 
speed for infinitesimal amplitude waves along the x direction 
in rubber as a function of X x and X y . In rubber, sound speeds 
vary strongly with the degree of extension due to nonlinear 
strain-extension relation and large deformations. In particular, 
when X x ^ X y , the sound speed depends on the direction. We 
use a time-of-flight method because standard ultrasonic tech- 
niques for measuring sound speeds are difficult to implement. 
The sheet is first stretched to an extension state X y > X x . Two 
record needles (marked by '+' in Fig. [3a)) are placed in con- 
tact with the sheet about 5 cm apart along the x direction. A 
bar is attached across the y direction of the sheet about 5 cm 
from the first needle. The bar is attached to a speaker and 
an in-plane, x-directed perturbation is applied to the sheet by 
pulsing the speaker. The signals from the record needles are 
recorded with a digitizing oscilloscope, and the longitudinal 
speed in the x direction is calculated from the time lag be- 
tween the separate signals and the known distance between the 
needles. Note that the speed is that of a linear wave moving on 
the stretched material; i.e., we are measuring the propagation 
speed of small distortions around the stretched state. 



We have been unable to measure the shear wave speed di- 
rectly due to the thinness of the material, which makes it diffi- 
cult to excite shear waves. Instead we measure the appropriate 
force-extension curve and calculated the velocity from it. As a 
proof-of-principle, we have calculated in this manner the lon- 
gitudinal wave speed. We measure the force vs. displacement 
curve for the configurations depicted in Fig. |3jb). The mate- 
rial is stretched to the state (X x ,X y ), a long but narrow portion 
is gripped along the long edges, and one of the grips is oscil- 
lated sinusoidally at / = 25 Hz in the x direction while the 
other was fixed. A load cell and accelerometer attached to the 
driven clamp measures the applied force F and acceleration 
a. The displacement Sx = a/(27r/) 2 follows by integration. 

Since the measurement is in a strip narrow compared to its 
length, departures of the strain in the y direction from the 
stretched state e yy is zero almost everywhere. This state of 
strain corresponds to the state excited by a longitudinal wave 
in a thin sheet: only e xx and e zz are nonzero. Therefore, the 
force and displacement measurements in this configuration 
yield the modulus needed to calculate cl- Y — FW/ (SxLd) 
where W is the distance between the grips, L is the length 
of the grips, d is the thickness of the material in the stretched 
state, and the longitudinal wave speed along the x direction is 
c l = y/Y/p, where p = 944 kg/m 3 is the mass density of the 
sheet. 
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FIG. 3: Experimental configurations for determining sound speeds: 
(a) direct measurement of the longitudinal wave speed using a time- 
of-flight method; indirect determination of (b) the longitudinal wave 
speed and (c) the shear wave speed from a force-displacement mea- 
surement. In each setup, the sheet is first brought to the desired state 
of extension (A^, X v ), the fixed and movable clamps are attached to 
the sheet, and the measurements proceed. 

The ratio of the measured to calculated wave speed is plot- 
ted in Fig. |3 for X y = 3.2 and X x varying between 2 and 
3.2. For perfect agreement this ratio should be 1. The data 
points lie within experimental error of this value, validating 
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the procedure for calculating the wave speed from the force- 
displacement relationship. 

Shear wave speed. The configuration in Fig. 3(c) is used to 
measure the shear modulus, G = FW/ (4LtSx), from which 
we calculate the shear wave speed in the x direction, cs — 
J Gj p. A rubber sheet is stretched to X y > X x , and a long 
narrow rubber strip is gripped on opposite edges, as shown. 
A thin bar down is glued to the center line of the strip and 
oscillated in the y-direction. The result for G obtained from 
the force and acceleration measurements then yields c s . 
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FIG. 4: The ratio of measured to calculated longitudinal wave speed 
plotted as a function of \ x for fixed \ y = 3.2. 

Material frame. It is interesting to examine our data in 
the material (Lagrangean) frame rather than the lab (Eulerian) 
frame because, as we show below, in the material frame the 
shear wave speeds are completely isotropic. In the material 
frame, length is measured in units of a coordinate system that 
deforms with the material: for biaxial deformation (A^, X y ), 
the relation between velocity in the material frame c'" 1 ' and 
the stretched (or lab) frame is c^™) = (c4 , ) = 

{ C X /Ajc, Cy 

Aw)- 

Results. In Figure |5] we compare crack speeds with wave 
speeds as a function of for X y = 3.2. The data indicate 
cracks travel 10-20% faster than the shear wave speed, but 
slower than the longitudinal sound speed. We adopt here the 
descriptor "intersonic" used in studies of shear loaded cracks, 
for cracks that travel at speeds between the longitudinal and 
shear wave speeds. 

Discussion There are some natural questions to ask con- 
cerning choices made in Fig. [5] (1) When rubber is stretched 
the wave speed is anisotropic. Why did we choose just to plot 
the wave speeds in the x direction? (2) Large amplitude de- 
formations typically travel at speeds different from those of 
infinitesimal amplitude waves. Why can we use small ampli- 
tude speeds as a basis for comparison? 

The answers to these question stem from the fact that 
our rubber samples are described accurately by the Mooney- 
Rivlin theory, which we now summarize. Any contin- 



FIG. 5: Comparison of measured longitudinal and shear wave speeds 
as a function of X x at fixed \ y = 3.2 with Mooney-Rivlin theory 
(dashed lines), (a) x(v) and y(») direction longitudinal wave speeds 
in the material frame, (b) Longitudinal wave speed (V), crack speed 
(•), and shear wave speed (A) in the material frame, (c) Longitudinal 
wave speed (V), crack speed (•), and shear wave speed (A) in the lab 
frame. 

uum theory of deformation that depends upon relative mo- 
tion of neighboring mass points can only involve parti- 
cle displacements x(X) through the strain tensor E a p = 
(l/2)[^ (dxj/dXadxj/dXp)— 8 a p\. Herex(X) describes 
the location in space of a mass point originally at X. A two- 
dimensional isotropic theory can only depend upon the strain 
tensor through the two invariant quantities I\ = E xx + E yy 
and I2 = E xx E yy — E xy . Under biaxial strain En = 
(A? - l)/2. In the Mooney-Rivlin Theory OH [3, the 
strain energy density of a thin sheet of rubber is given by 

W = a(h + bl 2 ) (1) 

In the material or Lagrangean frame a calculation of the lon- 
gitudinal and shear wave speeds yields: 

Pc x 7 )2 = a[l+ b -(Xl-l)} (2) 

pc ( v 7 )2 = ^ + \^l-l)\ (3) 
pc x ^=a(l- b -) (4) 
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To recover the speed in the lab or Eulerian frame c( s ) = 

(d s \cy S ^) — (X x Cx , XyC^). From these equations we see 
that the longitudinal and shear wave speeds in the x direction 
are independent of X x and that the y direction longitudinal 
wave speed increases monotonically with increasing X x , both 
of which we observe. More importantly, ci'"" 1 is independent 
of direction. With the value of a = 501 (m/s) 2 and b = 0.106, 
all our sound speed data are well characterized by these equa- 
tions, as is shown by the dashed lines in Fig. [5] 

Furthermore, in Mooney-Rivlin materials, speeds of in- 
finitesimal longitudinal and shear waves are exactly equal to 
the finite amplitude plane wave speeds. Therefore, although 
we are measuring small amplitude wave speeds, these speeds 
bound the speed of any finite amplitude plane wave deforma- 
tion. Thus, it is appropriate to compare the wave speeds we 
measured to crack speeds. 

Lastly, we note that when a moving object travels in an elas- 
tic medium faster than the wave speed, it creates a shock wave 
in the form of a Mach cone. The wedge-like crack opening 
typical of cracks in rubber (shown in Fig.|2j is strikingly simi- 
lar to the Mach cone, suggesting that rubber cracks are shock- 
like, exceeding some response speed of the medium. Shock 
waves are known to form when finite tensile pulses travel in 
pre-stretched rubber strips 1 17]. 

Conclusion. In earlier work we found that cracks in rub- 
ber follow oscillating paths when the horizontal extension ex- 
ceeds a critical value. We still have no good explanation for 
this dynamical transition. However, knowing that the frac- 
tures travel above the shear wave speed will be important in 
resolving this question. 
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